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ABSTRACT 



We present new high-resolution spectro-polarimetric Can /18542 observations of umbral flashes in sunspots. At nearly 0/18, and 
spanning about one hour of continuous observation, this is the most detailed dataset published thus far. Our study involves both LTE 
and non-LTE inversions (but includes also a weak field analysis as a sanity check) to quantify temperatures, mass flows and the full 
magnetic field vector geometry. We confirm earlier reports that UFs have very fine structure with hot and cool material intermixed 
at sub-arcsecond scales. The shock front is roughly 1,000 K hotter than the surrounding material. We do not observe significant 
fluctuations of the field in the umbra. In the penumbra, however, the passage of the running penumbral waves alter the magnetic field 
strength by some 200 G (peak-to-peak amplitude) but it does not change the field orientation (at least not significantly within our 
sensitivity of a few degrees). We find a trend of decreasing high-frequency modulation power for more inclined magnetic fields for 
the line-of-sight velocity and magnetic field strength. In the outer penumbra we find an absence of high frequency power while there 
is increasingly more power at high frequencies towards the umbra. 
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1. Introduction 

Umbral flashes (UFs) are sudden brightenings observed in the 
chromospheric core of the Can lines. They tend to appear with 
a periodicity of roughly 3 minutes (Beckers & Tallant 1969; 
Wittmann 1969) and it is now reasonably well established that 
they are the chromospheric counterpart of the photospheric os- 
cillation (Kneer et al. 1981; Thomas 1984). The waves, prop- 
agating from the underlying photosphere, steepen and become 
non-linear as the atmospheric density drops by several orders 
of magnitude. Eventually, the waves develop into shocks which 
are observable as blue-shifted emission reversals in the Ca n line 
cores. A semi-empirical study by Felipe et al. (2011) allowed to 
reproduce multi-line sunspot observations using a 3D non-linear 
MHD code (Felipe et al. 2010). A more theoretical study by Bard 
& Carlsson (2010) analyzed which modes contribute to the for- 
mation of UFs and investigated the role of the coronal radiation 
in the strength of these shock events. 

The development of spectro-polarimetry has paved the way 
for new advances in our understanding of UFs, mostly aided by 
the characteristic nearly anti- symmetric shape of the Stokes V 
profile. The first polarimetric observations of UFs (Socas- 
Navarro et al. 2000a,b) already revealed the occurrence of 
"anomalous" Stokes V profiles during the flash events. The 
anomalous profiles could only be explained as the superposi- 
tion of two regular profiles of opposite polarity and considerably 
Doppler-shifted with respect to each other. The obvious conclu- 
sion then is that, while UFs are a large-scale phenomenon, form- 
ing emission patches that are several arc-seconds across, they 
harbor fine-scale structure. Therefore, even inside a single pixel, 
one has a mixture in the horizontal direction where a certain fill- 
ing factor is occupied by a "quiet" component (no reversal and 



zero or slightly downflowing velocity) while the rest is occupied 
by the shock-wave (line core emission reversal and strong up- 
flows). 

Additional indirect evidence supporting the two-component 
scenario described above was presented by Centeno et al. (2005) 
using observations of the Hei 10830 A multiplet acquired with 
the TIP instrument (Collados et al. 1999). However, such fine 
structure inside the umbral flash patches was not directly seen 
until years later, with the analysis of Socas-Navarro et al. (2009) 
using high-resolution images of the Ca n H line core from the 
Hinode satellite (Kosugi et al. 2007; Tsuneta et al. 2008). With 
the increased spatial resolution offered today by some ground- 
based telescopes, it has been possible to confirm this finding by 
other authors (Bharti et al. 2013; Henriques 2013). 

Umbral flashes seem to be connected to the phenomenon 
known as running penumbral waves, an oscillatory pattern orig- 
inally observed in Ha images (but now also in Hei, as shown 
by Bloomfield et al. 2007, and Can in this work). The wave pat- 
tern is seen to move horizontally across the penumbra, running 
away from the umbra in the radial direction (Tsiropoula et al. 
2000; Rouppe van der Voort et al. 2003). The relationship be- 
tween UFs and penumbral waves is not only suggested by vi- 
sual appearance. The first spectro-polarimetric observations of 
Lopez Ariste et al. (2001) showed evidence that the shock bends 
radially outwards following the field lines as it propagates up- 
wards. More recently, the inversions of infrared Stokes spectra 
carried out by Bloomfield et al. (2007) indicate that the same 
wave propagation that gives rise to the umbral flashes is at work 
in the penumbra and produces the running penumbral waves in 
a scenario where the magnetic field is more inclined. 



de la Cruz Rodriguez et al: Physical Properties of a Sunspot Chromosphere with Umbral Flashes 



The paper is arranged as follows: the observations and data 
processing are described in Sect. 2. A detail description of the 
umbral flashes in our observations and a simple analysis of the 
magnetic field strength, using the weak-field approximation are 
provided in Sect. 3. In Sect. 4 we describe the non-LTE inversion 
code Nicole and our results are presented in Sect. 5 and Sect. 6. 
Finally, we summarize our conclusions in Sect. 7. 



2. Observations and data processing 

The data used for this work were acquired at the Swedish 1-m 
Solar Telescope (SST, Scharmer et al. 2003) using the CRisp 
Imaging SpectroPolarimeter (CRISP, Scharmer 2006; Scharmer 
et al. 2008). CRISP is a Fabry-Perot interferometer (FPI) con- 
sisting of two etalons mounted in tandem. With this kind of in- 
strumentation one is able to record quasi-monochromatic images 
at a given wavelength of choice. A line profile may then be con- 
structed by scanning the corresponding spectral range, obtaining 
images at each line position. 

The observations were obtained on 201 1 May 4 at 08:30 UT. 
The target was a sunspot (AR 11204) at heliocentric (x,y) = 
(-336", 332"), observing angle 30° (ji = 0.87). 

The dataset includes full-Stokes observations in Can /18542 
and non-polarimetric images in H i /16563 (Ha). The Can /18542 
line was observed at 14 line positions with a sampling of 75 mA 
between A A = ±300 mA (an effective integration time of 
900 ms), and points at AA = ±400, ±860, +3100 mA from line 
center (with an effective integration time of 210 ms). Images are 
acquired concatenating four polarization states. Each polariza- 
tion state is a linear combination of the for Stokes parameters 
(/, Q, U, V). The noise, relative to the continuum intensity, is 
cr(Q, U, V) = (0.09, 1.0, 1.3) x 10" 3 . The Ha line is observed 
at AA = +0, +1300 mA from line center. 

Fig. 1 illustrates the observed line positions at /18542, which 
are marked on a quiet- Sun spatial average from the surround- 
ings of the sunspot. Ideally, one point in the continuum would 
be desirable but the presence of secondary transmission lobes in 
the CRISP transmission profile restricts the observations to the 
central part of the prefilter (see the bottom panel in Fig. 1). 

The data processing was done following a similar ap- 
proach as presented by Schnerr et al. (2011) but with some 
differences in the flat-fielding process, further discussed in 
Appendix A. High spatial resolution and precise alignment 
between the sequentially recorded images for different tun- 
ings was achieved with the Multi-Object-Multi-Frame-Blind- 
Deconvolution (MOMFBD, van Noort et al. 2005) technique for 
image restoration. 

In certain areas, where the scale of the seeing motions is 
smaller than the assumed size of the iso-planatic patch in the 
MOMFBD code, residual distortions can be present within the 
images of a same scan, or even the different polarization states 
at a given wavelength. Reducing the patch size during the re- 
construction can help to minimize this effect, but in some cases 
it does not remove it entirely. Henriques (2012, 2013) proposed 
to include extra wide-band images (one per wavelength and po- 
larization state) that do not contribute to the wavefront- sensing 
in the MOMFBD process, but contain the same residual distor- 
tions as the narrow-band images. The distortions are measured 
by comparing those extra images with the anchor wide-band im- 
age and each narrow-band image is de- stretched accordingly. 

The polarimetric calibration of the instrument was per- 
formed in 2D, resulting in a different Mueller matrix for each 
pixel (van Noort & Rouppe van der Voort 2008). The telescope- 
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Fig. 1. The Can /18542 observations were acquired at the line 
positions indicated on the line profile in the top panel. The ex- 
posure time is 900 ms (circles) and 210 ms (crosses). The bot- 
tom panel illustrates the CRISP instrumental profile at 854.0 nm 
(solid gray) in logarithmic scale. A solar-atlas profile multiplied 
with the prefilter (dashed-dotted line) is indicated in blue. 



induced polarization was removed using the SST model devel- 
oped by Selbing (2010) with the parameters fitted at /18542. The 
telescope model is further described in Appendix B below. 

The time series has been compensated for telescope rota- 
tion, alignment over time and residual rubber-sheet deformations 
(Shine et al. 1994) using the wide-band images as reference. 

In order to normalize the data to a spatially-averaged quiet- 
Sun continuum, we computed the intensity ratio between [i = 
0.87 and \i - 1.0 at AA - +3100 mA from line center, using 
synthetic profiles, computed in non-LTE, from a 3D hydrody- 
namical simulation (de la Cruz Rodriguez et al. 2011). This ra- 
tio was then applied to the intensity at the same wavelength of 
the McMath-Pierce Fourier Transform Spectrometer atlas (here- 
after FTS atlas, Brault & Neckel 1987), which has been pre- 
viously convolved with the CRISP transmission profile. In this 
manner we impose that the spatially-averaged quiet-Sun profile 
at AA = 3100 mA (pi = 0.87), must have an intensity of 0.792 
relative to the continuum at pi - 1 . 

Fig. 2 summarizes our observations. In the photosphere (see 
the wide-band image in the top-right panel), photospheric granu- 
lation surrounds an irregular sunspot. It has developed penumbra 
in approximately 3/4 of its circumference and a light-bridge di- 
vides the umbra in two halves (umbral cores). There is also a 
small pore on the top-left corner of the image. 

The Ha chromospheric image (top-left) is dominated by 
very dynamic fibrils which extend outwards, almost radially 
from the umbra-penumbra boundary. Some fibrils originate in 
the fore-mentioned pore and are oriented towards the sunspot, 
suggesting some magnetic connection between the two. 

The lowermost four panels correspond to the four Stokes pa- 
rameters at -150 mA from line center at /18542. The Stokes / 
image is remarkably similar to that in the core of Ha. 
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Fig. 2. CRISP filtergrams at the core of the Ha line (top left). All other images are in Can /18542. A wide-band image showing the 
photosphere is displayed in the top right panel. The lowermost four panels show the chromosphere at AA = -150 mA from line core 
in Stokes / (middle left), Q/I (middle right), U/I (bottom left) and V/I (bottom right). The arrows on the wide-band image indicate 
the direction to the solar North and to disk center, respectively. The largest white box indicated in the wide-band image highlights 
the subfield that has been used for the non-LTE inversions, whereas the smaller square indicates a subfield used to study running 
penumbral waves. 
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The Stokes Q and U signals, produced by magnetic fields 
perpendicular to the our line of sight, are mostly observed on the 
limb side of the sunspot, preferentially along elongated fibrils 
that extend in the superpenumbra (Bray & Loughhead 1974). 
The Stokes V image (associated with the longitudinal magnetic 
field component) shows very strong dynamics, with frequent 
sign changes due to emission reversals. 

The time series shows that UFs and running penumbral 
waves are present in the sunspot on both sides of the light-bridge. 



3. Umbral flashes at 07 18 

The /18542 data presented here are, to our best knowledge, 
the highest-resolution spectroscopic observations of UFs. Fig. 3 
shows the detailed evolution of a subfield (marked in Fig. 2 with 
the white square in the lower right corner) for eight time steps 
of our series. In Stokes / (first and second row), signatures of 
the UF are visible from At = 32 s to At = 80 s, and the sunspot 
seems to be in a quiescent state everywhere else. During the UF, 
our images reveal dynamic fine structure within the umbra that 
resemble penumbral micro-jets (Katsukawa et al. 2007), simi- 
larly reported by Socas-Navarro et al. (2009). 

In Stokes V, the imprint of the UF is dramatic. During quies- 
cence, the Stokes V images show negative polarity everywhere 
in the umbra. As the flash starts at At = 32, the presence of po- 
larity changes is tightly correlated with brightenings in Stokes /. 

Our observations also show the imprint of micro-jets above 
the penumbra in Stokes V. They are seen very conspicuously 
when the time series is played as a movie. 

The Stokes Q and U signals are too weak inside the umbra. 
However, the fibrils that form the superpenumbra on the limb 
side of the spot, show strong signal in Stokes Q and U within 
A A = ±300 mA. Running penumbral waves are visible in the 
four Stokes parameters. In Sect. 6 we analyze them in more de- 
tail. 



3.1. The weak-field approximation for Zeeman-induced 
polarization 

The weak-field approximation provides a relatively simple way 
of computing the magnetic field vector from full-Stokes obser- 
vations (see Landi Degl'Innocenti & Landolfi 2004, and refer- 
ences therein). It works particularly well in those cases where 
the Doppler width of the observed line (AA&) is larger than the 
Zeeman splitting produced by the magnetic field (AA B ). Under 
these conditions, the observed Stokes Q, U and V profiles can 
be reproduced as a function of the magnetic field strength and 
the first derivative of the Stokes / profile. 

The Ca n /18542 line has a relatively low Lande factor (g = 
1.10) and it is broader than most photospheric lines. Therefore, 
the weak-field approximation can be safely used even when the 
magnetic field is relatively strong. The inequality AAd > AAb 
provides a convenient means to calculate the regime of applica- 
bility for the approximation (Asensio Ramos 2011): 
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where m e and eo are the electron charge and mass, Aq is the 
central wavelength of the line, c is the speed of light, k is the 
Boltzmann constant, M is the mass of the atomic element, T is 
the temperature and v m { c the microturbulence. 
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Fig. 4. A time-series at a single pixel within the sunspot is shown 
in the top three-panels, in Stokes / (top left) and V (top middle). 
The top-right panel corresponds to a synthetic Stokes V pro- 
file computed using the weak-field approximation. The resulting 
longitudinal magnetic field strenght is plotted in the second row. 
For clarity, the lowermost panels illustrate the Stokes / and V 
profiles at a quiescent state (black) and during the UF (blue). 



Assuming a chromospheric scenario with T = 4500 K, 
^mic = 3 km s" 1 and calcium atoms, the upper limit is approxi- 
mately B < 2500 G. Even though this value is lower than usual 
photospheric values in sunspots (B « 3000 G), the magnetic 
field strength is usually weaker in the chromosphere. 

One limitation of the weak-field approximation is that it also 
assumes a constant magnetic field. While most photospheric 
lines are sensitive to a relatively thin layer of the atmosphere, 
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Fig. 3. A time-series showing a subfield of the FOV in Stokes / (the two uppermost rows) and Stokes V/I (two lowermost rows) at 
A/l = -150 mA. The UF is clearly visible in the central part of the series (At = 48 - 80 s). The minor tickmark separation is Off 5. 
An online movie is available (http://folk.uio.no/jaime/movies/fig3_delacruzetal2013.mov). 



where the magnetic field is nearly constant, the /18542 line is 
sensitive to a range of approximately 1500 km (e.g. Cauzzi et al. 
2008) and it is less justified to assume a constant magnetic over 
that height range. Another possible approach would be to use 
the bisector method (Rayrole 1967; Landi Degl'Innocenti & 
Landolfi 2004), which allows one to retrieve some depth infor- 
mation but it is more sensitive to noise. 



In this work we use the weak-field approximation to infer an 
approximation to the line-of- sight (l.o.s.) magnetic field strength 



at chromospheric heights. To keep a simple approach, we use 
only A/l - ±225 mA from line center. 

Fig. 4 summarizes our results for one pixel in the sunspot 
as a function of time. The top row illustrates t - A diagrams for 
Stokes / and V. The rightmost panel shows the Stokes V profiles 
that result from our fits using the weak-field approximation. 

The figure shows nicely the Stokes V polarity reversals dur- 
ing the UF, which coincide with emission in the Stokes / profile. 
In fact, the shapes of the observed Stokes V profiles are tightly 
reproduced by the first derivative of Stokes / and the polarity re- 
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versals are consistent with the changes in the slope of Stokes / 
that are induced by the UF. 

The top-middle panel in Fig. 4 illustrates the inferred longi- 
tudinal magnetic field component during the time series. Despite 
some fluctuations with time, the magnetic field strength remains 
rather constant showing that it is possible to reproduce UF pro- 
files with a constant magnetic field. We have indicated all the UF 
phases using short vertical lines at the bottom of the panel. 

For clarity, we have selected one quiescent profile and one 
UF profile, marked in the top row with black and blue horizontal 
lines, respectively. Both profiles are plotted in the lowermost two 
panels. The corresponding time steps are marked in the upper- 
middle panel, showing that the field strength is very similar dur- 
ing the quiescent and flash phases (at least within the limitations 
of the weak-field approximation which provides us with a rough 
average of the magnetic field). 

These results confirm that the strong variability of the /18542 
line is likely to be produced by changes in the temperature and 
velocity of the plasma and, to a much lesser extent, by changes 
in the magnetic field. Note that the field is probably perturbed by 
the shock waves during the UF, but it is not a main contributor 
to that variability observed in Stokes / or V. 

4. The inversion code Nicole 

To derive atmospheric parameters from our spectropolarimetric 
observations, we use the non-LTE inversion code Nicole (Socas- 
Navarro et al. 2000c, Socas-Navarro et al. in prep). Nicole solves 
the non-LTE problem using preconditioning as explained by 
Socas-Navarro & Trujillo Bueno (1997), assuming that each 
pixel harbors a plane-parallel atmosphere. All Zeeman sub- 
levels originating from a given atomic level are assumed to be 
equally populated, discarding any quantum interference between 
them, as proposed by Trujillo Bueno & Landi Degl'Innocenti 
(1996). We use a Can model atom consisting of five bound 
level plus a continuum (see Leenaarts et al. 2009) and the colli- 
sional broadening is computed according to Barklem & O'Mara 
(1998). 

The velocity-free approximation is used to compute the atom 
population-densities. By neglecting the velocity-field, only half 
of the profile needs to be computed and fewer quadrature angles 
are used to compute the mean intensity at each depth-point. This 
approximation is justified by the very broad shape of the Can 
lines and the fact that the atomic level populations are dominated 
by the radiation field at the line core. Once the atomic popula- 
tions are converged, the emerging full- Stokes vector is computed 
including the velocity field using a quadratic DELO-Bezier for- 
mal solution of the polarized radiative transfer equation (de la 
Cruz Rodriguez & Piskunov 2013). 

Scattering polarization is not included in our calculations 
but this is not a problem because the magnetic field is suffi- 
ciently strong to produce Zeeman-dominated Stokes profiles. In 
the /18542 line, it is safe to assume Zeeman induced signal when 
the magnetic field is stronger than approximately 100 G (further 
details can be found in Manso Sainz & Trujillo Bueno 2010; 
Carlinetal. 2013). 

Nicole has been previously used in LTE studies (Socas- 
Navarro 2011; Scharmer et al. 2012) and validated in non-LTE 
using a chromospheric 3D MHD numerical simulation (de la 
Cruz Rodriguez et al. 2012). In that work, the synthetic pro- 
files are degraded with a realistic CRISP transmission profile. 
In order to properly convolve with the instrumental profile, the 
synthetic spectra are calculated first in a fine grid of regularly 
spaced wavelength points. The inversion code works internally 



Table 1. Summary of node points used during each cycle of the 
inversion. Note that Bi . s is oriented along the z-axis. 
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with such finely sampled profiles except when it comes to com- 
paring with observations, at which point only the observed line 
positions are used to drive the inversion. 

The inversion process is initialized with a smoothed version 
of the HSRA atmosphere (Gingerich et al. 1971). Node points 
are located equidistantly along the depth-scale, where correc- 
tions to the model are applied. The correction in between nodes 
is determined as a cubic Bezier-spline interpolation (see Auer 
2003) of the node values. 

To improve convergence properties, the inversion is run in 
two cycles as proposed by Ruiz Cobo & del Toro Iniesta (1992). 
In the first cycle, fewer degrees of freedom (nodes) are pre- 
scribed. Once the solution cannot be further improved, a new 
cycle is started with an (slightly) increased number of nodes. The 
degrees of freedom used for each cycle, summarized in Table 1, 
have been obtained after some trial-and-error experiments with a 
few sample profiles. Note that we have used fewer nodes than in 
de la Cruz Rodriguez et al. (2012) because our observations are 
slightly undersampled around line core, and very few points are 
present in the photospheric wings of the line. Therefore, some 
depth information is lost with respect to what could be obtained 
with more observables. 

The equation of state defines a unique relationship among 
temperature, electron pressure and gas pressure. If two of them 
are known, then the third one can be obtained from it. During 
the inversion, however, only the temperature is set for any guess 
model. Therefore, hydrostatic equilibrium is used to compute the 
gas pressure. 

A recent study by van Noort (2012) proposes to use the in- 
strumental point- spread-function (PSF) to couple the inversions 
spatially. His method improves the convergence rate of inver- 
sions, producing very smooth results and it allows to resolve 
details at the diffraction limit of the telescope. However, it re- 
lies on precise knowledge of the instrumental PSF, which is not 
always the case in ground-based observations, where the seeing- 
related PSF needs to be taken into account. Here we use the 
traditional 1.5D approach, where each pixels in inverted inde- 
pendently from the rest, and therefore less consistency and more 
noise is expected in the 2D maps of inferred quantities. 



5. Atmospheric properties during umbral flashes 

To characterize the atmospheric state during UFs, we have in- 
verted two time steps from the series, corresponding to At = 16 s 
and At = 64 s in Fig. 3. To accelerate the computations, only 
every second pixel in both spatial directions was inverted. We 
are not interested in the very small-scale (close to the telescope 
diffraction limit) details of the umbra. Therefore, this loss of res- 
olution (now 0724) does not influence our conclusions signif- 
icantly whereas computation speed is decreased by a factor of 
four. 
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Fig. 5. Top row: Observed (dashed line) and inverted (solid line) 
profiles from one pixel within the sunspot during the UF phase 
(blue) and quiescence (black). Middle row: Density map illus- 
trating the inferred temperature stratification during UF (left) 
and quiescence (right) for a patch inside the umbra. A total of 
280 models is shown. Bottom row: Same as middle row, for the 
l.o. s velocity stratification. Negative velocities correspond to up- 
flows. 



Fits to sample profiles are shown in Fig. 5. Observed (dashed 
line) and synthetic (solid line) profiles are plotted for a quiescent 
phase (black line) and during a UF (blue line) in Stokes / (left) 
and V (right). Stokes Q and U profiles are too noisy to obtain 
meaningful fits in the center of the spot. The inversion code is 
able to reproduce to a large extent the spectral features present 
in the observations, including the emission reversal seen at line 
center during the UF. The middle and bottom rows show the in- 
ferred quantities as a function of depth, within a l'/75 x 2'/ 35 
patch around the point marked with the cross in Fig. 6, which 
are discussed below along with 2D maps of the inferred quanti- 
ties. 

In interpreting inversion results, it is important to understand 
that we cannot "see" a discrete height. If we take an atmosphere 
and make a perturbation (e.g., increase the temperature) in a slab 
that is much narrower than the photon mean free path, then the 
perturbed slab would be completely transparent and the radia- 
tion would travel through it unchanged. Thus, the spectral pro- 
files emerging from this perturbed atmosphere would be indis- 



tinguishable from the original unperturbed one. What this means 
that it is impossible to determine observationally the properties 
of a very thin layer. We are sensitive only to the average proper- 
ties of the atmosphere on length scales that are long enough to 
leave an imprint on the outcoming radiation. This is not a limi- 
tation of the inversion technique or the instrumentation. It is an 
intrinsic limitation imposed by the integral form of the radia- 
tive transfer equation. How severe this limitation is in practice 
depends on a number of parameters, the most important being 
the number of spectral lines observed. Obviously, having several 
lines probing the same atmospheric region increases the depth 
resolution of the inversions. 

Because of the limitation in height resolution discussed 
above, it is not always sensible to interpret inversion results at 
a specific height. Instead, one needs to average in the vertical 
direction over a slab that is thick enough to have a measurable 
effect on the observed profiles. Therefore, the 2D maps presented 
in this section have been averaged over a depth range of 0.5 dex, 
centered at the log T500 value indicated in each of the panels. 

The temperature stratification within the umbra changes dra- 
matically from the quiescence to the UF phase. The results from 
our inversions are illustrated using 2D maps at several optical 
depths in Fig. 6. Note that a constant r slice usually maps a cor- 
rugated height surface. 

In the middle photosphere, at logrsoo = -0.5, there is not 
a significant difference between quiescence and UF (lowermost 
panels of Fig. 6). We can only see those intensity changes pro- 
duced by solar evolution. Some umbral dots are visible within 
the umbra, and part of the light bridge is included in the upper- 
left corner of the field of view (FOV). 

In the upper photosphere (middle panels), the light bridge 
and surrounding penumbra become smoother and most of the 
fine structure of umbral dots has disappeared. 

Signatures of the UF start to appear within the sunspot at 
approximately logrsoo = -3.0, where the temperature starts to 
increase relative to the quiescent phase. This temperature ex- 
cess is very conspicuous at log T500 ~ -4.5. The lower-left panel 
in Fig. 5 shows that the photospheric temperature is essentially 
identical during both phases, starting to diverge approximately 
at logTsoo = -3. The distributions of models in this small sub- 
field shows that the maximum temperature difference between 
both phases, about 1500 K, occurs at logTsoo = _ 4. 

The 2D temperature maps in the chromosphere (top row in 
Fig. 6), show a hot canopy above the umbra during the UF that is 
not present during quiescence. This canopy coincides with those 
areas showing enhanced line-core emission. 

When shocks are present, we have a very sharp discon- 
tinuity of the physical parameters in the vertical direction. 
Discontinuities are blurred out by the radiative transfer and, as 
discussed above, observations do not allow to distinguish a sharp 
change from a smoother transition, (de la Cruz Rodriguez et al. 
2012) demonstrated how in similar situations, and having only 
the /18542 line, the inversion retrieves a somewhat smoothed out 
version of the shock. 

Fig. 7 shows the l.o.s velocity component at logTsoo = -4.2 
and log T500 = -2 from top to bottom respectively. In the chro- 
mosphere (top row), the inversion code reconstructs upflowing 
velocities up to 5 km s" 1 in the flash (left panel). During the 
quiescence, the umbra shows an homogeneous downflow patch 
(~1 km s" 1 ), similarly to the scenario reported by Socas-Navarro 
etal. (2001). 

Another interesting result is found on the light bridge. The 
time series shows that micro-jets have apparent up and down mo- 
tions in the intensity images. Our maps show that those motions 
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Fig. 6. Inferred temperature from a non-LTE inversion. 
Temperature slices from the inversion are shown in two columns 
at quiescent (right) and flash (left) phases at several optical-depth 
logTsoo = -4.2, -2, -0.5 respectively from top to bottom. Note 
that we have performed a depth average of 0.5 dex, centered at 
those heights. The minor tick-mark separation is Of! 5. The white 
cross marks the center of the patch of which the temperature and 
velocity stratifications are shown in Fig. 5 



are in counter-phase with the UF within the umbra, suggesting 
that the flash has an effect on the motions of the light bridge, as 
well. 

One limitation of our observing setup is that we have only 
one pseudo-continuum point in the far wing of the line. This 
point allows us to probe the lower photospheric temperature but, 
lacking a symmetric point on the other side of the line, the ve- 
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Fig. 7. l.o. s. velocity component from a non-LTE inversion. 
Slices at constant optical depth are shown in two columns at 
quiescent (right) and flash (left) phases at logTsoo = -4.2,-2 
respectively from top to bottom. Negative velocities correspond 
to upflows. The minor tick-mark separation is Of! 5. 



locities at that height are poorly determined. This explains why 
the temperature map is smooth in the photosphere whereas the 
corresponding velocity map is much noisier. 

The longitudinal component of the magnetic-field recon- 
structed by the inversion is noisy. We have compared the results 
from inversions with 1 node and 2 nodes respectively in Bi ong . 
The resulting maps are almost identical when a depth average 
between logTsoo = _ 4.5 and logTsoo = _ 3.5 is considered, al- 
though we notice that the inversion noise increases when the de- 
grees of freedom is higher. 

Further testing with the inversion weights (different for each 
of the Stokes parameters) showed that the reduced wavelength 
coverage in our observations allows the code to over-fit the noisy 
Stokes V profiles, which are usually perfectly reproduced by the 
fitted spectra. 

Increasing the degrees of freedom makes this situation 
worse, therefore we have used the inversion including only 1 
node, illustrated in Fig. 8 and Fig. 9. The former illustrates the 
resulting B\ ong across the FOV, assuming a constant magnetic 
field as a function of depth. As expected, the magnetic field 
strength is stronger in the center of the umbra, and becomes 
weaker outwards, rather typical configuration for this kind of 
sunspots. Some fibril-like features are visible in the lower-left 
corner of the FOV, where the superpenumbra starts. 

Even though we do not find significant changes in the con- 
figuration of the magnetic-field (Bi ong ) during the flash, small 
scale variations of approximately 200-300 G are present within 
the umbra. Those variations would be compatible with a noise 
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Fig. 8. Longitudinal magnetic field component from a non-LTE 
inversion with one node in Bi ong , during the UF phase (left) and 
quiescence (right). The minor tick-mark separation is 07 5. 
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Fig. 9. A 2D probability map illustrating the relation between 
the longitudinal magnetic field component during the UF and 
quiescence, over the inverted field of view. 



of 10" 2 in Stokes V (a rough estimate computed with the weak- 
field approximation, the FTS atlas profile and a magnetic field 
of 2000 G), however, our observations reach 10" 3 in the range 
A A = ±300 mA, suggesting that the source of these variations 
may not be the noise. 

Additionally, small scale features are present over the light- 
bridge (top-left corner), indicating a less smooth and more con- 
fined magnetic field configuration, suggesting a connection with 
the chromospheric jets that here are also present over the light- 
bridge. 

Fig. 9 illustrates a probability diagram comparing Bi ong (XJF) 
vs. B\ ong (quiescence) over the inverted FOV. This figure is com- 
patible with an scenario where the magnetic field does not 
change much during the UF, as we could somewhat conclude 
from our simple analysis using the weak-field approximation. 

The lack of signal in Stokes Q and U inside the umbra does 
not allow for a relevant interpretation of the horizontal compo- 
nents of the magnetic field, however, we include 1 node in B x 
and in B y respectively so the code can compute a Zeeman split- 
ting compatible with the observations where the magnetic field 
strength is close to 2500 G (see Sect. 3). 



6. Running penumbral waves 

de la Cruz Rodriguez et al. (2012) justify and demonstrate the 
suitability of LTE inversions to derive the l.o.s. velocity and the 
magnetic field vector from chromospheric observations, disre- 
garding the retrieved "temperature" merely as a fitting param- 
eter parameter with no actual physical meaning. The idea is to 
allow the inversion code to adjust the temperature stratification 
in such a way that the resulting LTE source function is the same 
that would be obtained in non-LTE with the correct temperature. 
Normally, in Ca n inversions this results in an unrealistically cold 
atmosphere lacking a transition region. The advantage of this ap- 
proach is that one can then analyze longer time- series or larger 
FOVs as the computing time per inversion is significantly re- 
duced. 

Note that Beck et al. (2013) also perform LTE inversions to 
retrieve atmospheric parameters in the upper photosphere using 
the extended wings of the Can /13968 (H line), but those authors 
restricted themselves to photospheric quantities. 

We have used LTE inversions to study running penumbral 
waves in our data. With a similar setup to the one described 
in Sec. 4, we inverted a small subfield over 193 snapshots. The 
FOV has been selected to include a region with strong Stokes Q 
and U signal (thus improving our sensitivity to the transverse 
field component). The FOV is indicated in Fig. 2 (see the small 
white box on the WB image) . This 2'.'5 x 5" patch has been con- 
veniently selected such that the running penumbral waves prop- 
agate parallel to the longest dimension. 

The resulting magnetic field vector has been converted into 
the solar reference frame to compute the local inclination of 
the field (similarly done in Scharmer et al. 2012). We solve the 
1 80 deg ambiguity by selecting, at each pixel, the direction that 
makes the field point radially outwards in the local frame. In 
Fig. 10 and 11 we average the plotted quantities in the direc- 
tion perpendicular to the wave propagation (the small dimension 
of the box) with the aim of increasing the signal-to-noise ra- 
tio. The l.o.s. velocity and longitudinal magnetic field are depth- 
dependent. For those quantities we plot depth averages from 
logr 50 o = -3.5 to logr 50 o = -4.5. 

The top panels of Fig. 10 show the l.o.s velocity (left), mag- 
netic field strength (middle) and the magnetic field inclination 
(right) as a function of time. Waves are visible in these diagrams 
as a modulation pattern of alternating bright and dark features. In 
the penumbra (left part of the panels), these features are slanted 
which is indicative of wave propagation. These are the running 
penumbral waves which are so prominently visible in the Stokes 
/ and V movies. In the umbra (x > 5) the modulation pattern is 
more horizontal and is related to UFs. The slanted features are 
clearly visible in the l.o.s. velocity and, with somewhat lower 
contrast, in the longitudinal field as well. The slope of the fea- 
tures indicate a propagation speed of about 13 km s" 1 , which is 
higher but consistent with the measurements of Rouppe van der 
Voort et al. (2003) who analyzed running penumbral waves in 
Can H and K images. 

According to the results shown in the middle panel, the mag- 
netic field strength is perturbed by the running penumbral wave, 
showing peak-to-peak variations up to 200 G. For clarity, the 
time series along the white dotted line is shown at the bottom 
panel. 

The rightmost panel in Fig. 10 illustrates the evolution of the 
magnetic field inclination as a function of time. The magnetic 
field is almost vertical close to the umbra and it rapidly becomes 
more horizontal outwards. Some modulation may be discerned 
as a function of time, but it is very weak and the fore-mentioned 
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Fig. 10. x-t diagrams for the l.o.s velocity (left), the magnetic 
field strength (middle) and the magnetic field inclination (right) 
are shown in the top row. The umbra is on the right side of the 
panels, the outer penumbra on the left side (x = 0). The de- 
tailed magnetic field strength modulation for one spatial location 
(at x = 3.65, indicated with the white dotted line in the middle 
panel) is plotted in the bottom panel. The magnetic field inclina- 
tion has been smoothed with a 3 pixels full- width-half-maximum 
Gaussian kernel. 



diagonal slope in the pattern is not clear. We then conclude that 
the field inclination is not perturbed significantly within our de- 
tection threshold (a few degrees) as the wave propagates. This 
result suggests that the shock front has an effect on the field 
strength but not on its orientation. 

The FFT power spectra along each of the columns of the 
l.o.s velocities and magnetic field strength (see Fig. 10) are pre- 
sented in Fig. 11. We find that in both cases, higher frequencies 
are more dominant towards the umbra (right side of the panels) 
where the magnetic field is more vertical (as indicated by the 
average magnetic field inclination derived from averaging over 
columns of the right panel of Fig. 10 and marked on the top axis). 
Outwards into the penumbra (lower values of x), most power is 
concentrated in lower frequency modes, producing a slope in the 
power spectra diagram. 
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Fig. 11. FFT power spectra along each column for the l.o.s ve- 
locities and magnetic field strength presented in Fig. 10. The axis 
labeling at the top indicates the average magnetic field inclina- 
tion obtained through averaging over the right panel of Fig. 10. 
More inclined field is found in the outer penumbra (x = 0), more 
vertical field in the umbra (right side of the panels). The solid 
horizontal line indicates the 3 minutes mode, whereas the hori- 
zontal dashed line indicates the 5 minutes mode. 



7. Conclusions 

In this paper we present a novel analysis of state-of-the-art spec- 
tropolarimetric observations of a sunspot chromosphere. In par- 
ticular, we have focused on the study of the physical parameters 
of the chromosphere during UFs. 

A simple analysis of the Stokes V reversals that are observed 
during UFs show that, even with a nearly constant magnetic 
field, the profiles can be successfully reproduced. Despite the 
limitations of the weak-field approximation, it is possible to con- 
nect the variability of Stokes V with the changes of the intensity 
profile (Stokes /). 

The resulting model from our non-LTE inversions suggest a 
temperature excess up to 1000 K in the shocked chromosphere, 
relative to the quiescence phase. In the photosphere, the temper- 
ature stratification seems unaffected by the UF phenomenon. 

To reproduce the observed profiles, a steep velocity gradient 
is reconstructed at logrsoo ~ -3. Above, the velocity remains 
almost constant around vi . s = 4.5 km s" 1 . These velocities are 
similar to those from the simulation by Bard & Carlsson (2010) 
at approximately 1000 km height. 

In the photosphere, the line is very broad and our ob- 
servations have limited coverage of photospheric wavelengths. 
Therefore, the inversion is rather insensitive to large photo- 
spheric velocities, which rather are adjusted to reproduce the 
fore-mentioned velocity gradient. 

Our observations do not allow for a detailed interpretation of 
the horizontal component of the magnetic field inside the um- 
bra: the umbral Stokes Q and U signals are too weak to allow 
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for robust results. Therefore we only comment on the longitudi- 
nal component. The overall magnetic field configuration remains 
rather constant inside the umbra during the flash, but we find 
subarcsecond fluctuations that could be produced by the flashes. 

The situation is slightly different in the penumbra, though, 
perhaps because the field is more inclined (resulting in a stronger 
component perpendicular to the l.o.s.) or because it is weaker 
(or a combination of both reasons). We inverted a 2D subfield 
along the superpenumbra where running penumbral waves were 
clearly observed. The magnetic field vector inferred by the inver- 
sions is almost vertical at the boundary of the umbra, becoming 
more horizontal outwards. Further analysis of the entire time se- 
ries have allowed us to measure the imprint of these waves in 
the magnetic field strength, which seems to be modulated with 
a peak-to-peak amplitude some 200 G. Our results, however, do 
not show a similar fluctuation of the magnetic field inclination, 
which suggests that the wavefront is able to change the field 
strength but not its geometry. 

We find a trend of decreasing high-frequency modulation 
power for more inclined magnetic fields for the l.o.s. veloc- 
ity and magnetic field strength. In the outer penumbra we find 
an absence of high frequency power while there is increasingly 
more power at high frequencies towards the umbra. This obser- 
vation is interesting in the context of wave propagation along in- 
clined magnetic fields. Low frequency waves (/?-modes), that are 
normally evanescent in the low density chromosphere, can leak 
into the chromosphere along inclined magnetic fields (Suematsu 
1990; De Pontieu et al. 2004). This mechanism has been identi- 
fied as the driving mechanism of dynamic fibrils (Hansteen et al. 
2006; De Pontieu et al. 2007; Heggland et al. 2007). This mech- 
anism might well be at play for UFs and running penumbral 
waves too and our analysis demonstrates that the kind of obser- 
vations we present here are well suited to study this in more de- 
tail. We leave the extension of the inversions to a wider FOV and 
inclusion of a wider range of inclination angles that are needed 
to be more conclusive on this matter to future work. 
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Appendix A: The flat-field correction 

Especial care must be taken with the flat-fielding procedure. 
At the SST, the flat-field data are acquired while the telescope 
moves in circles around solar disk center, accumulating « 1000 
images per wavelength and polarization state. Effectively, the 
data represent a temporal and spatial average of quiet- Sun data 
at solar disk center. Therefore, it is reasonable to assume that the 
observed line profile at each pixel should be almost identical, ex- 
cept for the dirt, wavelength dependent fringes and the prefilter 
transmission, that can change across the FOV. 

The telecentric configuration used in CRISP is optimal to 
achieve high spatial-resolution observations, although the instru- 
mental transmission profile is not constant across the field-of- 
view (FOV). Pixel-to-pixel wavelength shifts (cavity errors) are 
present due to imperfections in the surface of the etalons that 
cannot be infinitely flat. Additionally, imperfections in the coat- 
ing of the etalons produce reflectivity variations across the FOV 
that change the width of the transmission profile (reflectivity er- 
rors). 



These effects become obvious when a spectral line is present 
at the observed wavelength, producing intensity fluctuations 
across the FOV that correlate tightly with the cavity errors and, 
less obviously, with variations in the reflectivity. 

While these instrumental effects are permanently anchored 
to the same pixels on the camera, solar features are moved 
around by differential seeing motions, however, the MOMFBD 
code is only able to process data where the intensity at any given 
pixel, is strictly changing due to atmospheric distortion. To over- 
come this problem Schnerr et al. (2011) proposed a strategy to 
characterize CRISP and prepare flat-field images that allow to 
reconstruct the data. 

Their idea is founded on the assumption that each pixel ob- 
served the same quiet-Sun profile during the flat-field data ac- 
quisition. To remove the imprint of the quiet-Sun profile from 
each pixel, an iterative self-consistent scheme is used to separate 
the different contributions to the intensity profile at each pixel, 
which we summarize in the following lines. To reproduce the 
observed profile on a pixel-to-pixel basis, they assume that the 
intensity at each pixel can be expressed as, 

hbsW = c • (I QS (A) * T(e ch , e rh )) • (1 + c 3 * A), (A.l) 

where I b s is the observed profile at a given pixel, Iqs is the esti- 
mate of the spatially-averaged quiet-Sun profile, T is a theoreti- 
cal CRISP transmission profile, cq is a multiplicative factor that 
scales the entire profile and C3 accounts for shifts of the prefilter. 
The * stands for a convolution operator. 

We simplify their strategy, by assuming constant reflectivi- 
ties across the FOV. Neglecting the effect of reflectivity errors 
works particularly well in broad-lines like the /18542, given that 
the transmission profile is relatively narrow compared to the 
Doppler width of the line. It also reduces the time required to 
compute the fits, from hours to minutes on a single CPU, as no 
convolutions are required any longer. We also allowed to include 
two more terms to the prefilter polynomial, given that at this 
wavelength there are fringes that change as a function of wave- 
length, introducing some distortion in the profiles. 

Fig. A.l (lower panel) illustrates a density map, where the 
profiles from all pixels have been divided by Co • (1 + c$A + C4A 2 + 
C5A 3 ) and shifted to the correct wavelength. The best estimate 
of the quiet- sun profile is plotted with a blue line. A tight cor- 
relation means that our model is able to reproduce the observed 
profiles, despite the simplifications applied to the original model 
by Schnerr et al. (2011). A cavity-error map is also shown in 
Fig. A.l 

A. 1. The CCD backscatter problem 

The back-illuminated Sarnoff CAM1M100 CCD cameras used 
in the CRISP instrument allow to acquire 35 frames-per- second, 
freezing atmospheric motions during each acquisition. These 
cameras work particularly well in visible light, but their effi- 
ciency decreases steeply in the infrared. 

Above approximately 7000 A, the CCD becomes semi- 
transparent and the images show a circuit pattern that is not 
present at other wavelengths, de la Cruz Rodriguez (2010) ex- 
plains a way to overcome this issue and properly apply the flat- 
field correction to the data. Furthermore, the images show a 
circuit-like pattern that cannot be removed by traditional dark- 
field and flat-field corrections, as illustrated in Fig. A.2 (top 
panel). 
Examination of pinhole data and flat- fielded data shows that: 
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Fig. A.l. Tbp: Map of the FPI cavity errors in units of wave- 
length shift. Bottom: Density plot of all pixels in the summed 
flat-field data (orange clouds) that have been shifted to the cor- 
rect wavelength and corrected for the different contributions in 
the model. The blue line is the average line profile described by 
a spline curve and closely resembles a disk center quiet Sun atlas 
profile. 

1. There is a diffuse additive stray-light contribution in the 
whole image. 

2. The stray-light contribution is much smaller in the circuit 
pattern and the gain appears to be enhanced. 

These problems can be explained by a semi-transparent CCD 
with a diffusive medium behind it, combined with an electronic 
circuit located right behind the CCD that is partially reflecting 
and therefore also less transparent to the scattered light. Under 
normal conditions an image recorded with a CCD, I , can be 
described in terms of the dark-field (Z)), the gain factor (G/) and 
the real image (7 r ): 
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Fig. A.2. Top: subfield from an individual acquisition, where tra- 
ditional dark-current and flat-field corrections have been applied. 
An electronic circuit pattern from the CCD chip is visible in the 
background. Bottom: same subfield, backscatter corrected. 
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Fig. A.3. sketch showing a conceptual model of the camera. The 
dark arrows indicate the incoming light beam from the telescope 
and the gray arrows represent the back- scattered light that re- 
turns to the CCD. 



However, in the infrared we need a more complicated model: 



I = D + f(l- f)[(G b G f I r ) * P]G b + fGflr 



(A3) 



I =D + (Gf I r ). 



(A.2) 



where / represents the overall fraction of light absorbed by the 
CCD, G b is the gain for light illuminating the CCD from the back 
which should account for the electronic circuit pattern. In the fol- 
lowing, we refer to G b as backgain. P is a Point Spread Function 
(PSF) that describes the scattering properties of the dispersive 
screen. In the backscatter term we assume that (1 - f)GfI r Gb 
is transmitted to the diffusive medium where it is scattered. A 
fraction of the scattered light returns to the CCD, passing again 
through the circuit. The sketch in Fig. A. 3 shows a schematic 
representation of the structure of the camera and the path fol- 
lowed by the light beam. 

To obtain the real intensity I r , the PSF P, the back-gain G b 
and the front-gain G/ must be known. The transparency factor 
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is assumed to be smooth because the properties of the dispersive 
screen seem to be homogeneous across the field-of-view. This 
allows us to include / in the front and the back gain factors 



1000 



Gy = Vwg, 

G f = f G f> 
so Eq. A.3 becomes: 

I = D + Gy[(GyG f I r ) * P] + Gfl r 



(A.4) 



This problem is linear and invertible, but a direct inversion would 
be expensive, given the dimensions of the problem. A numerical 
approach can be used to iteratively solve the problem. We define 



J = (GyG f I r )*P 



(A.5) 



In the first iteration, we initialize / assuming that the smearing 
caused by P is so large that / can be approximated by the average 
value of the observed image multiplied with the back gain. The 
back gain Gy is assumed to be 1 for every pixel in the first itera- 
tion. Furthermore, we assume that the product Gfl r can be esti- 
mated from Eq. A.2 ignoring back scattering, i.e., Gfl r ~ I -D. 
These values are of the same order of magnitude as the final so- 
lution, and therefore correspond to a reasonable choice of initial- 
ization. 

Jo * (GyG f I r ) * (Gy(I - D)\ (A.6) 

where the initial guess of the PSF P represents an angular aver- 
age obtained from a pinhole image. The small diameter of the 
pinhole only allows to estimate accurately the central part of 
the PSF. The wings of our initial guess are extrapolated using 
a power-law. The estimate of Gfl r is then 



Gflf — Iq — D — Gy J 



(A.7) 



which can be used to compute a new estimate /. This procedure 
is iterated, until Gfl r and I are consistent. This new value of 
/ is used to improve our estimate of Gy , by applying the same 
procedure to images that contain parts being physically masked 
(I r = 0), as in Fig. A.4. In the masked parts where I r = we 
have, 

I o -D-GyJ = 9 (A.8) 



so the back-gain can be computed directly: 



Gy 



Io-D 
J 



(A.9) 



Thus, every pixel must have been covered by the mask at least 
once in a calibration image in order to allow the calculation of 
the back-gain. With the new Gy we can recompute a new esti- 
mate of Gfl r . 

We now need to specify a measure that describes how well 
the data are fitted by the estimate of P and Gy(P). Since both 
Gfl r and Gy are computed based on self-consistency, we itera- 
tively need to fit only the parameters of the PSF. We assume that 
the PSF is circular- symmetric and apply corrections to the PSF 
at node points placed along the radius. 

We use Brent's Method described by Press et al. (2002) to 
minimize our fitness function. This algorithm does not require 
the computation of derivatives with respect to the free param- 
eters of the problem. When the opaque bars block a region of 
the CCD, an estimate of the back gain can be calculated for a 
given PSF according to Eq. A.9. In our calibration data, the four 
bars of width L are displaced 0.5 L from one image to the next. 




200 400 600 800 1000 

X [pixels] 

Fig. A.4. Calibration data acquired to infer the backgain and 
the PSF describing the backscatter problem. The dark areas are 
physically masked on the focal plane of the telescope. The color 
scale is saturated to enhance the CCD chip pattern. 



This overlapping provides two different measurements of the 
back gain on each region of the CCD. However, as the location 
of the bars changes on each image, the scattered light contribu- 
tion is different for each of these measurements of the back gain. 
Our fitness function minimizes the difference between these two 
measurements of the back gain. 

Having thus obtained the PSF P and the back-gain Gy , we 
obtain Gf by recording conventional flats and assuming I r is a 
constant in order to obtain Gf from Eq. A.4. 

Appendix B: The SST polarization model at /1854 nm 

The optical path of the SST, from the lm lens down to the cam- 
eras, contains optical elements that can polarize light. Therefore, 
polarimetric observations at the SST must be calibrated. 

At the SST this calibration is split in two parts: the optical 
table and the turret. The former is calibrated daily, given that 
the observer can easily place a linear polarizer and a retarder to 
generate known polarization states, and analyze how the optical 
table changes the polarization (see van Noort & Rouppe van der 
Voort 2008, for further details). 

The result is a modulation matrix M tab i e that allows to con- 
vert the observed polarization states I\ c = (IkoJiciJiaJic?) 7 
into Stokes parameters / = (7, Q, U, V) T parameters, using a lin- 
ear combination: 



r = m. 



1 j 

table ' A lc> 



(B.l) 



where the prime symbol indicates that only the optical table has 
been calibrated. 

In principle, if the telescope was built with non-polarizing 
elements then (/', Q\ U\ V) = (/, g, U, V), but this is not the 
case. The turret has mirrors and a lens that change the polariza- 
tion properties of the incoming light. In fact, the polarizing prop- 
erties of the turret change with the pointing. Therefore, a final 
transformation (M te i) is needed to obtain telescope-polarization 
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free Stokes parameters: 



/ = M 



■tel 



•^-M-J 



(M-L-//c) 



(B.2) 



Selbing (2010) studied the polarizing properties of the SST 
and proposed a model to characterize the telescope induced po- 
larization at 630 nm. In this work, we use a similar model. Their 
work allows to compute M te / based on the azimuth and elevation 
parameters of the telescope pointing. 

Calibration images have been recorded during an entire day, 
using a 1-m polarizer mounted at the entrance lens of the SST. 
The polarizer rotates 360 deg in steps of 5 deg and several frames 
are acquired for each angle. These data are used to derive the 
parameters of the model proposed by Selbing (2010) at 854 nm. 
This model is already presented in de la Cruz Rodriguez (2010) 
and included here for completeness. 

Each polarizing optical element of the telescope is repre- 
sented with a Mueller matrix. L represents the entrance lens. 
The form of the lens matrix represents a composite of random 
retarders, so it demodulates the light without polarizing it. The 
reference for Q is aligned with the lm linear polarizer axis and 
the values of the matrix are measured in that frame: 



L(lo,h,fa>h>h) 



(I 

h h -h 

h h U 

l 2 -k (/o + h~l)J 



Mirrors are noted with M and have two free parameters. 
Assuming a wave that propagates along the z-axis and oscillates 
in the x - y plane, the parameters are the de- attenuation term 
R between the xy components of the electromagnetic wave and 
the phase retardance S produced by the mirror. This form of M 
assumes that Q is perpendicular to the plane of incidence: 



M(R,S) = 



2 



~R r — tzR 

■jR ~ + ~R 





-Vfl cos(^- Q 6) 
VRsin(^i5) 



-Vflsin^S) 

■y/Rcos(^6)) 



Finally, R(a) corresponds to a rotation to a new coordinate 

frame: 

(\ (^ 

cos(^or) sin(^ar) 

-sin(^or) cos(^ar) 

,0 1. 



R(<*) 



Selbing (2010) builds the telescope model using the Mueller ma- 
trix of each polarizing element, as a function of the azimuth (cp) 
and elevation (0) angles of the Sun at the time of the observation 
(Eq. B.3). We have included the conversion factor from degrees 
to radians in the matrices, thus all the angles are given in degrees: 



M te i (0, <p) =R f+ (c n ) ■ M s (c 9 , cio) • M f (c l9 c 8 > 

• Rf-(c n ) ■ Raz(<P) ' M az (c 5 ,c 6 ) • R el (0) 

• M e i(c 5 , c 6 ) • L(c , cuc 2 , c 3 , c 4 ), 



(B.3) 



where ci,2,...,n are the parameters of the model, computed us- 
ing a least-squares fit to the observed data. The fitted values are 
summarized in Table B.l. 

Note that the polarizer has a significant leak of unpolarized 
light at 854.2 nm. Using small samples of the sheets used to con- 
struct the 1-m polarizer, we have estimated the extinction ratio 
of the 1-m polarizer to be approximately 40%, so the parameters 
of the lens have been determined assuming that value. 

The linear polarization reference is defined by the first mirror 
after the lens, however this is not very useful in practice because 



Table B.l. Parameters of the telescope model obtained from cal- 
ibration data at 8542 nm. 



Parameter 


Fitted value 


Co 


+9.596161 


io- 1 


C\ 


-2.541791 


lO" 3 


Cl 


-8.337375 


IO" 3 


C3 


+9.640901 


10" 1 


ca 


-1.875988 


IO" 2 


C5 


+9.164327 


10" 1 


ce 


+ 1.465364 


10 1 


Cl 


+ 1.017471 


10° 


c& 


+ 1.326977 


10° 


c 9 


+ 1.010035 


10° 


ClO 


-1.592178 


10° 


C\l 


+4.000000 


10° 



the turret introduces image rotation along the day. Instead, we 
use the solar north as a reference for positive Q by applying an 
extra rotation to M te i. 

The full-Stokes Can images in Fig. 2 have been calibrated 
using the telescope model described in this section. 
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